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In the past, measurement results of splice loss of optical fibers have 
corresponded poorly to existing theory, which assumes a uniform power 
distribution across the cone of radiation defined by the local numerical 
aperture. In this paper, a model is developed in which a Gaussian power 
distribution across the local numerical aperture is assumed. Trans- 
mission through a splice at each point on the transmitting core is found 
to depend on the ratio of receiving to transmitting numerical aperture 
at that point. Numerical integration of these "point" transmission 
functions over core areas of interest yields both splice loss and the 
additional loss that occurs in a long fiber following the splice. This 
model cannot be theoretically rigorous, since it is inconsistent with 
boundary conditions required by the laws of light propagation. How- 
ever, it has been found to predict splice loss under varying conditions 
with much greater accuracy than existing theory. The model has the 
further virtue of being able to calculate how variations in many intrinsic 
and extrinsic splice parameters combine to produce an overall splice 
loss. 

I. INTRODUCTION 

Calculations for the loss in an optical fiber splice, as a function of 
offset, tilt, diameter, or numerical aperture mismatch, have been re- 
ported by several authors. 1-3 These calculations all assumed a uniform 
power distribution across the cone of radiation defined by the local nu- 
merical aperture (NA). This is consistent with assuming equal mode 
excitation, equal mode attenuation, and no mode coupling. 4 An as- 
sumption concerning the power distribution is necessary to characterize 
all combinations of both intrinsic and extrinsic splice imperfections. 
While these assumptions allow easy calculations, correspondence with 
measurement data has been unacceptable. 3 Gloge 5 reported results based 
on a diffusion process from the uniform power distribution to the 
steady-state distribution calculated by Marcuse. 6 Correspondence with 
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measurement data is much improved for the case of small offsets; how- 
ever, approximations used in the theory make it difficult to characterize 
other splice imperfections or splice loss for large offsets. 

This paper presents a phenomenological model with an assumed power 
distribution selected solely due to resulting correspondence of calculated 
effects of splice imperfections with measurement data. The model is not 
intended for uses other than splice loss characterization, and since the 
model does not obey the laws of ray optics exactly, caution should be 
exercised in any other uses of it. 

II. DETAILS OF THE MODEL 

The power distribution across the cone of radiation defined by the 
numerical aperture at any point on a fiber core is assumed to be Gaussian 
in form. Figure 1 is a sketch of the assumed radially symmetric distri- 
bution across the cone of radiation from a given point on the fiber core. 
The model consists of solving for the transmission ratio in terms of 
transmitting and receiving numerical apertures, then integrating this 
ratio over the core areas of interest. 

Consider the steady-state power distribution across the cone defined 
by the local numerical aperture to be Gaussian in form, normalized to 
a value of 1 at r = 0. 



p( r ) = e-'W 
where a is proportional to the width of the Gaussian. 



(1) 




Fig. 1 — Gaussian power distribution. 
3168 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1978 



Numerical aperture has been defined by various investigators to be 
the 1/e 2 , 0.1 or 0.01 power level. This power level, po, determines the 
width of the Gaussian and will be left as a parameter. Therefore, at 

Po = e -roW (2) 

or p(r) = e {r2/r ° 2)ln p°. (3) 

III. NA MISMATCH 

The application of this simple model can best be illustrated by con- 
sidering the change in the power distribution as it propagates through 
a splice. The usual class of circularly symmetric index of refraction 
profiles 4 is used throughout this paper. In this example, transmitting 
and receiving fiber profiles are identical except for the value of »oi. 
Therefore, the NA as a function of transmitting fiber core radius (p) is 

NAi(p) « n i V2E[ [ 1 - (£) " 1 1/2 > (4) 

where Ai « (moi — re c )/noi is small 

«oi = refractive index at center of core 
n c = refractive index of cladding 
a — 2 for nearly parabolic profile fibers 
R = fiber core radius. 

The angle <j>i(p) (Fig. 1) is determined by NAi(p) to be 

0i(p) = sin" 1 [NAi(p)/m(p)] f (5) 

where 

ni(p) = n i[l-2A 1 ^) a | L/2 . (5a) 

The power distribution at a given value of p is assumed to be 

Pl ( r )= e (r2/ri2)lnpo 5 ( 6 ) 

where ri is proportional to tan 4>\ (Fig. 1). 

For the receiving fiber, NA2(p) will first be assumed to be less than the 
transmitting fiber NAi(p), so that 

p 2 (r) = e (r2/r 2 2 > ln P<> (7) 

is the relative power distribution that corresponds to the same point on 
the receiving fiber. Figure 2 shows these two functions for r<i < r\ and 
for Gaussian distributions which are truncated at r 2 and r\, respectively. 
Using these truncated distributions, we assume the power in region I is 
lost immediately at the splice, since power in this region lies totally 
outside the receiving NA. Splice transmission through this point is 
then 
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Fig. 2 — Point transmission distributions. 



t( P ) = 



i 



Pi(r)rdr 



o ^ 



(8) 



r) rdr 



t(p) = — — [ e Wn)*m po - i] (9) 

Po- 1 
for r 2 < ri. 

Since r\ = k tan <pi and r 2 - fe tan 2 where fe is a constant of propor- 
tionality, then using (4) and (5), 

r 2 tan (sin" 1 V2A 2 ) _. NA 2 .^ 

n tan (sin -1 V2Ai) NAi 
for some given p. 

For this particular example, NA2/NA1 is constant for every point (p) 
on the fiber core; therefore, eq. (9) gives the total transmission just after 
the splice for a truncated Gaussian distribution. 

For the case of nontruncated Gaussian distributions, (8) and (9) be- 
come 



t( P ) = 



J 2 p 1 (r)rdr+ f p 2 (r) rdr 
Jri 



I 



Pi(r) rdr 



£(p) = l+(^) 2 po-e (r 



2/ri) 2 ln po 



(ID 



(12) 
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for r 2 < ri. For small values of p , (12) and (9) yield similar results. The 
full Gaussian is used to simplify later computations. 

Since this example contains radially symmetric distributions, there 
is no dependence on 0. Equation (12) gives the transmission through the 
splice at a point for the steady-state power distribution assumed in (3). 
Again, for this example, (r2/Vi) 2 is constant across the fiber core; there- 
fore, (12) gives the total transmission just after the splice for a numerical 
aperture mismatch. 

IV. COMPARISON WITH NA MISMATCH DATA 

Equation (12) can be compared to measurement data for a splice with 
an NA mismatch (r 2 < r\, receiving fiber NA < transmitting fiber NA). 
This transmission coefficient represents an immediate loss at the splice. 
Figure 3 contains measurement data along with calculations for both 
a uniform power distribution and a Gaussian distribution. The param- 
eter, p , is set to 1/e 2 and 0.1. As shown in Fig. 3, sensitivity to the value 
selected for pn is significant only for large NA mismatches. 

Figure 3 contains data obtained using a HeNe and a GaAlAs laser 
source with a long input fiber (>500 m). Fibers selected for these mea- 
surements were well matched in O.D., core diameter, and a, but con- 
tained mismatches in NA. There were, however, slight differences in O.D., 
core diameter, and a on the order of a few percent, so that measured loss 
would be expected to differ somewhat from calculated values. 
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Fig. 3 — Loss due to numerical aperture mismatch. 
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Correspondence with measured data using either source is good 
compared to the uniform power distribution model which predicts too 
much loss. 

Before this model can be applied in general, the case of transmission 
from a smaller NA transmitting fiber to a larger NA receiving fiber must 
be considered. If, in Fig. 2, r^ represents the transmitting fiber and r\ 
the receiving fiber, then all the power contained in the cone of radiation 
defined by the transmitting NA is within the receiving NA. A unity 
transmission coefficient is assumed for this case. 

V. NEAR-FIELD POWER DISTRIBUTION 

The near-field power distribution as a function of core radius, p, can 
be calculated for the Gaussian model after a suitable weighting function 
is applied. Since a Gaussian distribution implies higher loss for higher 
order modes of propagation (the tails of the distribution), a weighting 
function is needed to reduce the amplitude of the Gaussian as a function 
of radius. This is required since only higher order modes of propagation 
are significant near the core-cladding interface. The weighting function 
assumed is the power distribution for the uniform power model. From 
(4), 

W"*[l-G) f ] W 

for the case of uniform power across the cone of radiation defined by the 
numerical aperture where Po is proportional to A, no and input power. 3 
With this function used as the amplitude at r = 0, then 

P(p)«P (l-jb) f (""eWr^PordrdO. (14) 

Using (4) to obtain ro, 

P(p)~PoU-P 2 ) 2 (15) 

where all constant terms have been combined in Po and R = 1. 

VI. COMPARISON WITH NEAR-FIELD POWER MEASUREMENTS 

Measurements of near-field power were made with a GaAlAs laser 
source after propagation through a long (>1 km) fiber wrapped under 
tension to simulate the effects of some microbending loss (~1 dB/km). 
These measurements were made using a 100X objective and a TV vidicon. 
camera. Figure 4 is a photograph of the camera output for a typical 
Western Electric fiber. An approximation to the curve was obtained by 
smoothing over the index dip and averaging power measurements for 
each side of the distribution. Calculations using (13) for the uniform 
distribution and (15) are plotted for comparison. The calculated 
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Fig. 4 — Near-field power distribution. 

steady-state power distribution using the Gaussian model with the 
(1 — p 2 ) .weighting function is in excellent agreement with the measured 
distribution. 

VI. ADDITIONAL LOSS IN THE FIBER AFTER A SPLICE 

A parabolic index fiber splice is known to cause additional loss in the 
fiber after the splice. 7 This additional loss depends on the differential 
mode attenuation and mode coupling characteristics of the receiving 
fiber and, for the fibers used in these experiments, is approximately equal 
to the loss at the splice for small offsets. 7 Therefore, any realistic de- 
scription of an optical fiber splice must include this effect. 
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Referring to Fig. 2, the power in region II is seen to be within the nu- 
merical aperture of the receiving fiber; however, this power is improperly 
distributed. The sharp decrease in power at the edge of the receiving 
numerical aperture, r<t, is physically impossible; however, this effect is 
significantly reduced compared to the uniform power model. As this 
distribution propagates down the fiber length (1), we assume that a new 
Gaussian steady-state distribution will be generated, as shown in 
Fig. 5. Some portion, c, of the excess energy contained in region II of 
Fig. 2 will be lost in the process of reestablishing steady-state conditions. 
If energy couples to adjacent modes with equal probability and if all 
modes are equally excited and attenuated, then 1/2 of the power in region 
II would be lost during redistribution. Since the excess energy in region 
II is skewed toward higher order modes which may be lossier, c would 
be expected to be greater than 1/2. The loss mechanism is probably 
transfer of some of this excess energy to higher order propagating or leaky 
modes. 8 

Referring to Fig. 2, the equivalent transmission loss due to power lost 
from region II is 



M(p) = 



j pi(r) rdr - j p 2 (r) rdr 
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(16) 
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Fig. 5 — Steady-state power distribution. 
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Evaluating this integral yields 

M(p) = -ceWrjptn po + c + c /^\ 2 .( Po - i). (17) 

To approximate c, we use the centroid for the solid of revolution for 
region II in Fig. 2. This places one-half the excess energy in the region 
< r < c and the remaining energy between c <r <r 2 . This radius is 
taken as the value for c. From Fig. 2, 

f C \piir) rdr - p 2 (r) rdr] = f * [pi(r) rdr - p 2 (r) rdr]. (18) 

Solving the integrals yields the following nonlinear equation: 

e (c/n)2[ n po _ /£S\ e (c/r 2 ) 2 ln po 

= -\l + e <r2/n)*i»po - p) 2 (po + 1)1. (19) 

This equation was solved for c/r 2 with < r 2 /ri < 1 for p = 0.1 and 
Po = l/e 2 , and found to be accurately approximated by a quadratic 
equation. For po = 0.1, 

c « 0.7994 - 0.08796 (— ) + 0.00846 (— ) (20) 

and for po = l/e 2 , 

c * 0.8041 - 0.0724 /-) + 0.00636 (— ) . (21) 

The total effect of the splice is then 

t toi (p) = t(p)-At(p). (22) 

itot(p) = (1 - c) [l + po 2 - «WW*» "»] + c0 2 , (23) 

where t{p) is the immediate point transmission coefficient at the splice 
given by (12) and At(p) is the reduction in transmission due to power lost 
in reestablishing a steady-state Gaussian distribution in a long fiber after 
the splice given by (17). Again, we call attention to the fact that the value 
of c probably depends on the differential mode attenuation and mode 
coupling characteristics of the receiving fiber. 

VIII. DIAMETER MISMATCH 

The necessary parts of the model have been developed so that the 
effect of diameter mismatch (Fig. 6) can now be considered. Let the 
transmitting NA function equal 

NAi(p) = no V2A \l - (j^) 2 j (24) 
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Fig. 6 — Core index of refraction profiles for the case of diameter mismatch. 

and the receiving NA function equal 

NA 2 (p) = n V2A Tl - (■£-) I L 2 ' 

where R\ and #2 are the transmitting and receiving core radii. For K = 
R2/R1 and R2 normalized to unity, 

2 i_ p 2 



(25) 



©-£ 



KV (26) 

This ratio is not constant with p (except for K = 1), so that an integration 
over the receiving core is necessary. Distributions remain radially sym- 
metric; therefore, no angle dependence is present. 



7\ot - 



C 1 tip) [I -K*fl*pdp 

J o [l-XV] 2 pdp 



(27) 



where t(p) equals eq. (12) for the effect at the splice, or by £ t ot(p) [eq. 
(23)], to include the additional loss in the fiber after the splice. This in- 
tegral must be solved by numerical techniques. 

IX. COMPARISON WITH DIAMETER MISMATCH DATA 

Figure 7 compares the calculations of short-length-diameter mismatch 
effects at the splice with measured results for various values of K. Fibers 
were drawn from the same preform to insure that a and NA mismatch 
were minimized. A HeNe laser source was used with approximately 1 
m of fiber after the splice. 

The Gaussian power distribution model shows good agreement with 
measurement data. Long-length-diameter mismatch data are not 
presently available. 
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Fig. 7 — Loss due to core diameter mismatch. 

X. TRANSVERSE OFFSET 

Perhaps the most important extrinsic splice parameter is transverse 
offset (axial displacement). The Gaussian model can be applied to the 
case of transverse offset; however, some computational difficulties are 
encountered. Referring to Fig. 8, the area of overlap is divided into re- 
gion I, where the receiving NA is greater than the transmitting NA, and 
region II, where the receiving NA is less than the transmitting NA. In 
region I, t(p) is unity, therefore with R = 1, 

J-»cos _1 d/2 /-»1 
J (1 - p 2 ) 2 pdpd<t> 

J d/2 

l'l = 



'd/2 cos <p 



J-» tt/2 •» 1 
I (1- P 2 ) 2 pdpd0 
Jo 



(28) 



In region II, 



ld/2 /» 1 



T„ = 



J-» cos 'd/2 /» 
Jd 



d/2 cos <j> 



t(q)[l - q 2 ] 2 pdpd<f> 



J» tt/2 r» 1 
| (l-q 2 ) 2 pdpd<t> 
o Jo 

where q = p 2 - 2pd cos + d 2 . 



(29) 
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Referring to Fig. 8, we note that the line separating regions I and II 
is a locus of equal NA; therefore, t(p) is unity on this line. On the curved 
boundary of region II, t(p) is zero; therefore, the derivative of t(p) is in- 
finite at the intersection of these lines. This point must be omitted and, 
since t{p) is steep in the vicinity of this point, a large number of incre- 
ments are required to evaluate eqs. (28) and (29) numerically. 

Ttot - Ti + T„. (30) 

If the loss at the splice is desired, (12) is used for t(p), and if the total 
loss including losses required to reestablish steady state is being calcu- 
lated, then (23) is used. 

XI. COMPARISON WITH TRANSVERSE OFFSET DATA 

Figure 9 is a comparison of transverse offset data and calculations 
(i) reported by Gloge, 5 (ii) using the Gaussian distribution assumption, 
and (Hi) using the uniform power assumption. Measurements were made 
with a GaAlAs laser source and an unbroken long fiber (2 km). After a 
reference level was established, the fiber was broken approximately in 
the center and the ends spliced to obtain the reference level again. Offset 
was introduced and the loss measured 1 km and 2 m after the splice as 
a function of transverse offset. Figure 9 shows the results for small offsets. 
Agreement among the Gloge model, the Gaussian model, and measure- 
ment data is good for loss at the splice. The Gaussian model also agrees 
well for the loss occurring in the fiber after the splice. 

Figure 10 shows results for large offsets. As compared with actual 
measurements, the Gaussian model understates the loss through a long 
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Fig. 8 — Regions of overlap for offset fiber cores. 
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Fig. 9 — Loss due to transverse offset. 

fiber after the splice for large offsets. The Gaussian model cannot, at 
present, account for this effect for large offsets. 

XII. CONCLUSIONS 

A point transmission model has been used to calculate splice loss 
due to NA mismatch, diameter mismatch, and transverse offset in 
parabolic-profile fiber splices. Near-field power distributions have also 
been calculated. Correspondence with measurement data is much im- 
proved as compared to calculations using the uniform power distribution 
model. 

The authors again emphasize that the Gaussian power assumption 
does not obey the laws of ray optics exactly and that this assumption was 
made primarily due to good correspondence of resulting calculations with 
measurement data and due to computational considerations. Future 
work will include the effects of a mismatch and the effects of long fibers 
after splices with combinations of intrinsic and extrinsic parameter 
mismatches. 
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Fig. 10 — Transmission vs offset for large offsets. 
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